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The experimental data on hadron yields and ratios in central lead-lead and gold-gold collisions 
at 158 AGeV/c (SPS) and y/s — 130 AGeV (RHIC), respectively, are analysed within a two-source 
statistical model of an ideal hadron gas. A comparison with the standard thermal model is given. 
The two sources, which can reach the chemical and thermal equilibrium separately and may have 
different temperatures, particle and strangeness densities, and other thermodynamic characteristics, 
represent the expanding system of colliding heavy ions, where the hot central fireball is embedded 
in a larger but cooler fireball. The volume of the central source increases with rising bombarding 
energy. Results of the two-source model fit to RHIC experimental data at midrapidity coincide with 
the results of the one-source thermal model fit, indicating the formation of an extended fireball, 
which is three times larger than the corresponding core at SPS. 

PACS numbers: 25.75.-q, 24.10.Pa, 25.75.Dw 
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I. INTRODUCTION 



Searching for the quark-gluon plasma (QGP) is one of 
the major objectives in the study of relativistic heavy ion 
collisions. The principal question is whether the strongly 
interacting nuclear, or rather parton, matter reaches the 
stage of chemical and thermal equilibrium. Although this 
idea was put forward by Fermi 50 years ago [1], up to now 
there is no unambiguous test to probe the degree of equi- 
libration in the system. One of the possible approaches 
is to study the equilibration process within microscopic 
models [2, 3, 4, 5, 6, 7]. The more traditional way is to fit 
macroscopic observables, yields and transverse spectra of 
particles, obtained in experiments to the statistical model 
(SM) of a fully equilibrated hadron gas. The simplic- 
ity of the SM has led to a very abundant literature (see 
[8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22] and 
references therein). As the experimental data became 
more precise, it has been understood that the ideal SM 
does not provide an adequate description of all hadron 
multiplicities [9]. Particularly, the yields of pions are 
usually underestimated while the abundances of strange 
particles are overpredicted. Besides, the calculated par- 
ticle number density and the total energy density at the 
stage of chemical freeze-out are found to be too large 
[15]. Therefore, some modifications to the SM have been 
proposed. These improvements consider the following ef- 



fects: (i) excluded volume effects that lead to a Van der 
Waals type equation of state (EOS) due to a non-zero 
"hard-core" radius of hadrons [13, 15, 16]; (ii) strangeness 
suppression, that enters into the distribution functions 
via the phenomenological factor 7s < 1 [11, 12]; (iii) 
chemical non-equilibrium of light quarks [18]. It is worth 
noting that the assumption of a single expanding source 
remains the basic ad hoc hypothesis of these models. The 
scenario of several fireballs has not been employed ex- 
cept a particular case in which all fireballs have the same 
chemical potential and the same temperature [17, 23] but 
move with different collective velocities along the beam 
direction. In this special case the particle ratios are, how- 
ever, not affected by superimposing a collective longitu- 
dinal kinetic energy and are, therefore, identical to the 
results of the single fireball scheme. Obviously, if the 
baryon density and/or the strangeness density are not 
the same everywhere in the reaction volume, the scenario 
with several independent sources cannot be reduced to 
the single source scenario. Our investigations of the two- 
source scenario have been inspired by the experimental 
observation of decreasing of antiproton to proton ratio 
with rising rapidity in the center-of-mass system [24] (the 
same is true for A/ A ratio also) which results to low net 
baryon densities in the midrapidity range of relativistic 
heavy-ion collisions at SPS energies (158 AGeV/c) [25]. 
In the two-source model the total reaction volume is di- 
vided into two regions, the inner source and the outer 
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source. Each source is assumed to be in the chemical and 
thermal equilibrium, respectively, and allowed to have 
different temperature, baryon and strangeness densities, 
etc. This two-source scenario is also in accordance with 
microscopic model calculations. Such microscopic cas- 
cade calculations show that the central zone between the 
remnants of the colliding nuclei is more baryon dilute 
compared to the baryon-rich zones of the nuclei residues 
[5, 6]. Thus, global equilibrium is not reached by the 
whole system, but local equilibrium in the central zone 
and in the peripheral region may still occur separately. 

The paper is organized as follows. Description of the 
two-source model is given in Sect. II. Sections III and 
IV present results of the model fit to the SPS and RHIC 
data, respectively. Comparison with the predictions of 
the standard thermal model is also given. Finally, con- 
clusions are drawn in Sect. V. 



II. TWO-SOURCE MODEL 

Firstly, we briefly sketch some of the basic principles 
of the statistical model of an ideal hadron gas. In the 
framework of the grand canonical ensemble, the macro- 
scopic characteristics of the system are derived via a set 
of distribution functions (we work with the system of 
units where c = Ti — ks = 1) 



f(p,m,i) 



exp 



\Jp 2 ~+rtf - VbB 1 



UsSi) IT 



il}- 1 , 



(1) 



where p and rrii are the momentum and the mass of the 
hadron species i, T is the temperature, ps and us are 
the baryon chemical potential and strangeness chemical 
potential, B% and S% are the baryon charge and strange 
charge of hadron i. The sign "+" in Eq. (1) stands for 
fermions and sign"—" stands for bosons. The particle 
number density n*, the energy density Si and the pressure 
P in the system read as 



n, 
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(3) 
(4) 



with gi being the degeneracy factor of hadron i and V 
the volume of the system. The entropy density s = S/V 
is determined via the Gibbs thermodynamic identity 



(e + P- HbPb - PsPs) T~ 



(5) 



Here ps and ps are the baryon density and strangeness 
density, respectively. 

In the present paper we develop a two-source statisti- 
cal model (TSM) of a hadron gas. The model divides the 



whole reaction zone into two regions: the outer region 
(source 1 or SI) and the inner region (source 2 or S2). 
Since these sources are assumed to be in local thermal 
and chemical equilibrium, they are allowed to possess dif- 
ferent temperatures, net baryon densities, chemical po- 
tentials, etc. Furthermore, the strangeness density is no 
longer kept zero everywhere as in the SM. Only the total 
number of strangeness is required to be zero because of 
strangeness conservation in strong interactions. 

The characteristics of a single fireball can be described 
by means of four independent parameters, such as vol- 
ume V, total energy E, net baryon density ps and net 
strangeness density ps- These parameters enable one to 
determine the fireball temperature T, the strangeness 
and baryon chemical potentials ps and ps- Since the 
total strangeness is zero in heavy-ion collisions, the num- 
ber of free parameters in the SM of an ideal hadron gas 
is reduced to three, namely, V, E (or e = E/V) and 
ps = Nb/V. In the two-source model the number of 
free parameters increases to seven. Although the net 
strangeness in each of the sources can be nonzero, they 
are linked via the the condition of total strangeness con- 
servation: 



N S i + N S2 = . 



(6) 



Our choice of the free parameters in the two sources is 
as follows: Besides temperature, volume, and net baryon 
density of the two sources, i.e. Ti, V\, ps\ and T 2 , V2, 
PB2, the strangeness density in source 1, psi, is consid- 
ered as a free parameter. 



III. HADRON PRODUCTION AT SPS 

The baryon yield and ratios of hadrons at midrapid- 
ity in central Pb+Pb collisions at 158 AGeV are listed 
in Table I together with the results of the TSM fit. 
All hadrons including their resonances and possible anti- 
particles with masses less than 2 GeV/c 2 are included in 
the fitting procedure. No additional constraints such as 
strangeness suppression or excluded volume are assumed 
except the feeding-back effect from resonance decay. For 
the sake of comparison, we check our model by fitting 
the same set of experimental data with the single-source 
model. The results of the fit are also listed in Table I. 
As in the TSM case, these results are obtained includ- 
ing the feeding-back effect but without further correc- 
tions. Compared to the ideal SM, the TSM improves the 
agreement with the experimental data. The thermody- 
namic quantities obtained from the two fits to experimen- 
tal data are shown in Table II. One can find the sizable 
differences in temperature and volume between the two 
sources. With 155 MeV the temperature in source 2 is 
significantly higher than that of 117 MeV in source 1. 
The volume of the hot source is about 0.7 Vo, where Vq 
is the volume of a lead nucleus, while the volume of the 
cooler source is four times larger. This fact indicates 
that the two-source object can be interpreted as a hot, 
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TABLE I: Baryon yield and hadron ratios at midrapidity for 
central lead-lead collisions at SPS energies and predictions of 
the single-source and two-source statistical models of an ideal 
hadron gas. 





Data 


SM 


TSM 


TSM-Data 0/ 
Data /0 


Ref. 


N B (net) 


372±10 


371.9 


372.1 





[26] 


K+/K- 


1.85± 0.1 


1.97 


1.87 


4 


[27] 


p/p 


0.07± 0.01 


0.069 


0.058 


-17 


[28] 




n 94Q+ n ni q 


9^1 


fl 947 


_l 




n/n 


0.383± 0.081 


0.427 


0.405 


6 


[29] 


A/A 
V /n° 
K° s /n~ 


0.128± 0.012 


0.130 


0.137 


7 


[30] 


0.081±0.013 


0.133 


0.108 


33 


[31] 


0.125±0.019 


0.121 


0.120 


-4 


[32] 


K° s /h- 


0.123±0.02 


0.102 


0.102 


-17 


[33] 


A/h~ 


0.077±0.011 


0.069 


0.064 


-17 


[33] 


Q/E 


0.219±0.045 


0.131 


0.104 


-53 


[33] 


S-/A 
X 2 /DOF 


O.llOiO.Ol 


0.156 


0.115 


4 


[33] 




46/9 


16/5 







TABLE III: The same as Table I but for the data set from 
Ref. [15]. 





Data 


SM 


TSM 


TSM-Data 0/ 
Data /0 


Ref. 


N B 


372±10 


362.9 


373.9 




f96l 


hr 


680±50 


606 


659 


-3.2 


f34l 


K° 

Ji s 


68±10 


61.6 


64.4 


-5.5 


[34] 


4> 


7.6±1.1 


13.4 


7.89 


4 


[35] 


p-p 


155±20 


125 


147 


-5 


[34] 


K+/R- 


1.80± 0.1 


1.99 


1.74 


-3 


[37] 


p/p 


0.07± 0.01 


0.065 


0.069 


-1.5 


[281 

L 




0.249± 0.019 0.220 0.229 


-8 


[29] 


n/n 

x 2 /dof 


0.383± 0.081 


0.411 


0.421 


10 


[29] 




41/6 


3.1/2 






T[MeV] 

V[fm s ] 

e[MeV/fm 3 ] 




157 

4160 

422 


114 (SI) 
155 (S2) 
17000 (SI) 
1920 (S2) 
60 (SI) 







298 (S2) 



TABLE II: Thermodynamic characteristics of two-source and 
single-source models. 





TSM (SI) 


TSM (S2) 


SM 


T[MeV] 
V[fm 3 ] 


117 


155 


158 


7250(3Vo) 


1705(0.7%) 


4203 


pslfm- 3 ] 


0.048 


0.014 


0.088 


Ps[fm-' A ] 


0.0018 


-0.0075 


0.0 


elMeV/fm 3 } 


104 


314 


436 


P[MeV/fm 3 ] 


14 


52 


67 


N B 


348 


55 


408 


Nb 





30 


36 


N B -Ng 


348 


25 


372 


N s 


34 


42 


117 


% 


21 


55 


117 


N s -Nj 


13 


-13 





p B [MeV] 


460 


45 


213 


ps[MeV] 


68 


17 


22 



relatively small core surrounded by a cooler and larger 
halo. The net baryon density in the halo is found to be 
14 times as large as that in the core. Hence, the major 
part of baryons is contained in the outer source, while the 
inner source contains almost all antibaryons. The total 
net baryon charge equals, however, the initial number of 
participants. 

The strangeness density is negative in S2 and pos- 
itive in SI. It means that the inner source contains 
more s-quarks than s-quarks. This finding is supported 
by microscopic model calculations [5]. From the mi- 
croscopic point of view the possible explanation of the 
phenomenon is as follows: Due to the strangeness con- 
servation in strong interactions strange and anti-strange 
particles must be produced in pairs. At SPS energies 
strangeness is mainly carried by kaons. Because of the 
small interaction cross section with hadrons, K + and K° 
are leaving the central reaction zone easier than strange 



particles which carry s quarks, e.g., A and K, thus lead- 
ing to a negative strangeness density in the midrapidity 
range. 

The energy density in S2 is about three times larger 
than that in SI. Such a low energy density in the 
outer source corresponds to the energy density at ther- 
mal freeze-out rather than at chemical freeze-out [5]. In 
other words, the solution for two sources can not be re- 
duced to the one-source picture even in the case where 
exclusive midrapidity data have been used. 

For further investigation we applied the two-source 
model to fit the set of experimental data used in Ref. [15], 
which are a combination of midrapidity and 47r-data. 
The data and the results of the TSM and SM fit are listed 
in Table III. Two important facts can be learned from 
this comparison. First of all, it is already known that 
the ideal SM (without strangeness suppression and ex- 
cluded volume effects) underpredicts the number of neg- 
atively charged hadrons h~ and overestimates the yield 
of mesons. The problem can be cured by increasing 
the volume V and decreasing the temperature T of the 
source [15]. However, the antibaryon to baryon ratios 
become then completely wrong. The TSM enables one 
to get both the correct multiplicities of h~ and <f) and the 
correct antibaryon/baryon ratios. Secondly, the temper- 
atures of both sources remain almost unchanged when we 
have shifted from the data set of Table I to the data set 
of Table III, and the temperature of the central fireball 
is very close to the temperature of the entire system in 
the single-source model. 

We have also fitted the hadron yields and ratios taken 
in the whole available rapidity range, see Table IV. 
Again, the thermodynamic characteristics in S2 are far 
from the averaged characteristics of the combined S1+S2 
system. For instance, the entropy per baryon, S/A = 
s/pb, equals 31 in source 1 but increases up to 373 in 
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TABLE IV: The same as Table III but for the hadron yields 
and ratios in the whole rapidity range. 





Data 


SM 


TSM 


TSM-Data 0/ 


Ref. 


N B 


372±10 


372.7 


373.5 


0.4 


[26] 


br 


680±50 


684 


673.4 


-1 


[34] 


K° 


68±10 


72.5 


68.8 


1 


[34] 


<t> 


7.6±1.1 


7.5 


7.6 





[35] 


p-p 


155±20 


138 


143 


-8 


[34] 




1.10± 0.05 


1.0 


1.0 


-10 


[36] 


K+/K- 


1.80± 0.1 


1.84 


1.80 





[37] 




0.125±0.019 


0.121 


0.118 


-5.5 


[32] 


X 2 /DOF 




1.3/5 


0.6/1 






T[MeV] 




126 


117 (SI) 












155 (S2) 






V[fm 3 ] 




17400 18500 (SI) 












1310 (S2) 






e[MeV/fm 3 ] 




90.6 


63 (SI) 












327 (S2) 







source 2, mainly due to the low net baryon density. The 
average entropy per baryon (Si + S2)/(A\ + A2) is 37 
which is very close to the value 37.8 obtained from the 
one-source SM fit. Results of all three fits favour the idea 
of the formation of a compact hot baryon-dilute central 
zone with the following averaged characteristics: temper- 
ature T = 157 ± 2 MeV, volume V = 0.6 ± 0.1 V , and 
baryon chemical potential ps = 31 ± 14 MeV. The tem- 
perature of the halo is much lower, i.e., Ts ± = 117 ± 3 
MeV. 



IV. HADRON PRODUCTION AT RHIC 

Experimental data on hadron yields and ratios in the 
midrapidity range in central gold-gold collisions at y/s = 
130 AGeV became available recently [38, 39, 40, 41, 42, 
43, 44, 45]. These data are listed in Table V together 
with the predictions of the SM and TSM. Surprisingly, 
now the two-source model divides the total volume in two 
parts with approximately equal sizes, temperatures, and 
other thermodynamic characteristics. The results of the 
SM fit and TSM fit are almost identical. It seems that 
the volume of the central fireball significantly increases by 
the transition from the SPS energies to the RHIC ones, 
and that hadrons detected in the midrapidity region are 
originated from a single thermalized source. Its volume 
is more than 5000 fm 3 , that is about three times larger 
than the corresponding core volume in Pb+Pb collisions 
at the SPS, and the temperature reaches 186 MeV. Such 
a high temperature value is close to the value T = 190 
MeV obtained in [42], while being 10 MeV higher than 
the temperature T = 175 MeV obtained in [22]. However, 
if the multiplicity of negatively charged hadrons, h~, is 
excluded from the data set, the temperature of the mod- 
elled system drops to 176 MeV in our calculations. This 
important question should be clarified in future studies. 



TABLE V: Hadron multiplicities and ratios at midrapidity 
for central gold-gold collisions at RHIC energy (y/s = 130 
AGeV) and predictions of the single-source and two-source 
statistical models of an ideal hadron gas. 





Data 


SM TSM 


TSM — Data 0/ 
Data 70 


Ref. 


N B 


343±11 


340.6 340.6 


-7 


[38] 


br 


2050±250 


2238 2239 


9 


[38] 


p/p 


0.60± 0.07 


0.57 0.57 


-4 


[39] 


p/ir~ 


0.08± 0.01 


0.087 0.087 


8 


[40] 


K-/K+ 


0.87± 0.08 


0.75 0.75 


-14 


[41] 


K'/tT- 


0.149±0.02 


0.153 0.153 


3 


[41] 


K* /hr 


0.060±0.017 0.036 0.036 


-40 


[42] 


K*°/h~ 


0.058±0.017 


0.030 0.030 


-48 


[42] 


A/A 


0.77± 0.07 


0.69 0.69 


-11 


[42] 




0.82± 0.08 


0.80 0.80 


-2 


r * ol 

[43] 


X 2 /DOF 




9.6/7 9.6/3 






T[MeV] 




185.7 185.9 (SI) 
185.5 (S2) 






V[fm 3 ] 




5157 2608 (SI) 
2560 (S2) 






e[MeV/fm 3 ] 




1297 1305 (SI) 
1282 (S2) 






fi B [MeV] 




52.6 53.3 (SI) 
51.8 (S2) 







We checked also that the incorporation of the excluded 
volume effects by assigning the hard-core radius r = 0.4 
fm to all particles leads to an enlargement of the total 
volume but does not affect the temperature of the fire- 
ball. 

Another interesting characteristics at the chemical 
freeze-out are the energy per particle and the entropy 
per baryon. In [17] the criterion E/N w 1 GeV was in- 
troduced for all particles, baryons and mesons, in a broad 
range of bombarding energies spanning from a few GeV 
per nucleon up to several hundred GeV per nucleon in 
the center-of-mass system. The predictions of the SM 
and the TSM for the midrapidity range are as follows: 
E/N = 1.1 GeV, and S/A = s/p B = 121. The last value 
is about 20% below the value s/pb = 150 predicted by 
the UrQMD calculations, but for Au+Au collisions at full 
RHIC energy, yfs = 200 AGeV [46]. It would be interest- 
ing to perform these microscopic model calculations also 
at Vs = 130 AGeV. 

V. DISCUSSION AND CONCLUSIONS 

In summary, the analysis of hadron multiplicities in 
central lead-lead and gold-gold collisions at 158 AGeV 
and y/s — 130 AGeV, respectively, is presented within a 
two-source statistical model of an ideal hadron gas. The 
TSM fit to the experimental data taken at midrapidity 
at RHIC practically coincides with the standard single- 
source fit. This result supports the idea of a formation 
of an extended hot fireball in the central zone of heavy- 
ion collisions at RHIC energies. The temperature of the 
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central fireball varies from 176 MeV to 185 MeV, depend- 
ing on the incorporation of the multiplicity of negatively 
charged hadrons, h~ , in the fitting data, whereas the ex- 
cluded volume effects seem not to affect the fireball tem- 
perature. Further going studies are necessary to answer 
this question. 

At SPS it is found that the properties of the system at 
chemical freeze-out can be well understood in terms of 
two sources, a central core and a surrounding halo, both 
being in local chemical and thermal equilibrium. Tem- 
peratures as well as baryon charge and strange charge 
of the two sources are different. The thermal character- 
istics of the central fireball, obtained from the TSM fit 
to hadron yields and ratios, depend only weakly on the 
considered rapidity interval of the data, i.e., midrapidity 
or the whole rapidity range. 

It is worth mentioning that strangeness seems to be in 
equilibrium in both sources which is reflected by 75 = 1 
in our calculations. This observation is in line with the 
fact that there is no need to introduce a strangeness sup- 
pression factor into the standard SM if one fits the parti- 
cle ratios from the midrapidity range of Pb+Pb collisions 
at SPS energies [16]. But the factor 75 < 1 arises if one 
intends to fit 47r-data [12, 15]. A possible explanation 
for this puzzle is a non-homogeneous distribution of the 
strange charge within the reaction volume. Therefore, 
the local, not global, equilibrium of strangeness can be 
reached separately in the central and in the outer part of 
the expanding fireball. 

Furthermore, in the standard SM the energy and 
hadron number densities are too high to treat the sys- 
tem at chemical freeze-out as a gas of point-like particles 
anymore. To resolve this problem the introduction of a 
repulsive hard-core potential for hadrons which leads to 
a Van der Waals type EOS might be important [15, 16]. 
In the case of two thermalized sources neither the energy 
density nor the hadron density is so large and the incor- 
poration of excluded volume effects (at least for the halo) 
becomes less important. 



Another interesting fact is that the temperature of 
source 2 at SPS energies is about 10 MeV lower than 
the temperature predicted by most of the single-source 
statistical models. In the SM such a low temperature 
can be obtained under the assumption of chemical non- 
equilibrium (overpopulation) of light u and d quarks [18]. 
In the latter case the system is not in a state with maxi- 
mum entropy, i.e., it is probably produced directly from 
an exploding quark-gluon plasma (QGP) and does not 
undergo further chemical equilibration in the hadronic 
phase [18]. However, the ratios of heavy strange baryons, 
such as fi/E or fi/E, are still underpredicted. The TSM 
allows for another interpretation of the data. It is quite 
likely that hot nuclear matter, which forms the core, 
spends enough time in the hadronic phase to reach a 
chemically equilibrated state. In this case the system 
does not remember about the past, e.g., the QGP stage 
(see also [47]). This circumstance complicates the de- 
tection of the plasma. It would be very interesting to 
investigate the forthcoming data on Au+Au collisions at 
RHIC energies (y/s = 130 and 200 AGeV) in the whole 
rapidity range within the TSM in order to check (i) the 
increase of the volume of the central fireball at the ex- 
penses of the halo; (ii) equilibration of strangeness in 
both sources; (iii) a possible change of the halo temper- 
ature. 
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